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Abstract 

The production of electroweak bosons (photons, W and Z particles) in PbPb and pp 
collisions at sqrt(s) = 2.76 TeV per interacting nucleon pair has been measured with 
the CMS detector at the LHC. Direct photon production is studied using samples of 
isolated photons. W and Z bosons are reconstructed through their leptonic decay 
into muons. Their production rate in PbPb data is studied as a function of the cen- 
trality of the collision and compared to that in pp interactions, once normalized by 
the number of binary nucleon-nucleon interactions. The results are also compared to 
next-to-leading-order perturbative QCD calculations. 
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1 Introduction 

The measurement of electroweak bosons (photon, W and Z) represents one of the basic tests of 
the Standard Model of elementary particles. In the context of heavy ion collisions these bosons 
are not expected to interact with the hot and dense strongly interacting medium created. In 
that case, their production rate and properties can be taken as a reference for other particles 
which result affected or provide normalisation for other physics processes. Thus, the direct 
comparison of the production cross section of such probes in pp and in nuclear collisions allows 
one to estimate possible modifications of nuclear parton density functions (PDF) relative to a 
simple incoherent superposition of nucleon PDF. 

Photons have already been studied in heavy ion collisions at several experiments, at the Rel- 
ativistic Heavy Ion Collider (RHIC) [1J and, more recently, also at the Large Hadron Collider 
(LHC) EJ|3l, but it has been the energy and luminosity conditions from the LHC that has made 
the detection of W and Z bosons possible. Studying these bosons in their leptonic decay modes 
turns out to be very useful, as leptons suffer negligible energy loss in the medium. Results 
on weak bosons production in PbPb collisions have already been published by ATLAS |U and 
CMS EHU collaborations. 

Several challenges are faced when studying the electroweak bosons in PbPb collisions. In the 
case of photon production, the main proccesses generating direct prompt photons (quark-gluon 
Compton scattering and quark-antiquark annihilation) are difficult to disentangle from the con- 
tribution of gammas from parton fragmentation processes, mainly 71° and rj decays, giving rise 
to a wide variety of produced photons. 

In the case of weak (W and Z) boson production, the main challenge lies in their low production 
cross section, compared to other processes. Although the leptonic decay modes are useful in 
identifying the boson production, they suppose an additional reduction factor (< 10%) of the 
yields. 

The CMS experiment [7] is well suited to detect muons and photons. In particular, for the 
photon detection and measurement, the electromagnetic calorimeter (ECAL) possesses a fine- 
grained (Atj x AO = 0.0174 x 0.0174 in < 1.44) hermetic disposition, arranged in a quasi- 
projective geometry, based on lead tungstate crystals which provide very narrow and precise 
electromagnetic showers. 

The silicon pixel and strips tracker allows tracking and transverse momentum particle mea- 
surement in the pseudorapidity range \r/\ < 2.5, with a good momentum resolution between 
2-3% for charged tracks of up to 100 GeV, thanks to the high magnetic field (3.8T) provided by 
the CMS superconducting solenoid. 

Muons are measured in the range \rjf\ < 2.4. Their triggering and identification takes place in 
the muon chambers system, embedded in the steel return yoke of the magnet, with an excellent 
performance, while the muon transverse momentum is measured in the silicon tracker. 

The centrality of the collisions is determined from the energy deposited in the forward steel/ quartz- 
fibre calorimeters. The measurements presented here are based on PbPb and pp data samples, 
both taken at y^s = 2.76 TeV per interacting nucleon pair, corresponding to an integrated lu- 
minosity of around 7 fib^ 1 and 231 nb _1 , respectively. Similar measurements have also been 
performed in larger samples of pp collisions at y 7 ^ =7 and 8 TeV, although they are not reported 
here. 
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2 Isolated Photons 



2 Isolated Photons 

Photons coming directly from the hard scattering of incident nucleons are usually highly ener- 
getic and produced isolated from other particles. This is not the case for photons from back- 
ground processes, mainly coming from decays of neutral mesons produced inside a jet, and 
thus, surrounded by significant hadronic activity from other parton fragments. Isolation pro- 
vides a good handle to suppress the fragmentation photon component in the data sample, 
while removing a very small contribution of direct photons. Still, the hard scattering is su- 
perimposed on top of other multiple parton-parton scatterings (underlying event) occurring 
simultaneously and which must be subtracted, prior to applying the isolation criteria. 

The contribution of the underlying event is subtracted, based on the mean value of the energy 
deposited per unit area along the r\ — <p strip containing the photon candidate in the ECAL. 
Once the mean energy is removed, the isolation requirements Ej < 5 GeV inside a cone 
defined by AR= 0.4, considering the energy deposited in the electromagnetic and hadronic 
calorimeters and Tracker) is applied. The resulting sample of isolated photons is dominated 
by quark-gluon Compton photons, with still a small proportion of photons from decays and 
fragmentation processes, which fulfill the isolation criteria. 

The shape of the shower deposited in the ECAL, given by the transverse width of the cluster 
(second moment of the electromagnetic energy cluster distribution around its mean tj position), 
is then exploited. Direct photons have a smaller cluster width, while background photons tend 
to have wider widths. A binned maximum likelihood fit of the data to two contributions is 
performed: one accounting for the signal photons (modelled with Pythia [8] embedded in PbPb 
data) and another one for the background, taking a template from the PbPb data. The number 
of isolated prompt photons is extracted from the fit. 

Figure [l] (left) shows the isolated photon yield per unit of photon transverse energy, Ex, de- 
rived from PbPb data and normalised to the number of binary collisions (N co u) by means of 
the nuclear factor ^a/^\ This quantity is presented as a function of the photon Ex and for 
three different regions in centrality ( [0-10%], [10-30%] and [30-100%]) as well as the whole 
minimum bias sample. The cross sections are extracted for the photon pseudorapidity region 
\t] 7 1 < 1.44, the ECAL barrel region. The corresponding cross section obtained in pp data, at 
the same centre-of-mass energy, is also included. The photon energy distributions from PbPb 
data, scaled by Taa, agree quite well with those from pp collisions and with NLO JETPHOX [0 
calculations. 

From this comparison, nuclear modification factors, Raa = t7 pbPb l (T')/ cr Pp(T') are derived, 
where c^pb (t) indicates the photon yield in PbPb data, divided by the normalising factor 
Taa- This is done for different centrality regions, as a function of photon Ex. Figure [l] (right) 
shows the [0-10 %] case, where the data points are compared to several model predictions for 
the nuclear PDF (EPS09 ED, nDS mi and HKN07 |T3|). The error bars account for the statisti- 
cal uncertainty and the coloured-filled rectangles for the systematic ones, which come mainly 
from the background template modelling and the photon energy scale, both in PbPb and in pp 
data. The Raa values are compatible with unity and, thus, confirm the scaling of the isolated 
photon production with N co q. The present uncertainties preclude any quantitative constraint 
on the (small) nuclear effects predicted by the different theoretical models. 

The dependence of the R^ factors with the number of participants in the collision, N par t/ for 
different photon Ex bins, can be seen in Figure |2j 

lr ^AA is the number of elementary binary nucleon-nucleon (NN) collisions, N co u, divided by the elementary NN 
cross section. 
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Figure 1: Left: Isolated photon spectra as a function of the photon Ej for several centrality 
regions in PbPb collisions (scaled by Taa) an d pp collisions, scaled by arbitrary factors (for 
better visualization). The results are compared to the NLO JETPHOX calculations. Right: Nu- 
clear modification factors, Raa, as a function of the photon Ej for the 0-10% most central PbPb 
collisions at \/s = 2.76 TeV . 
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Figure 2: Nuclear modification factors, Raa, as a function of PbPb centrality (given by the 
number of participant nucleons, N part ) for five different photon Ej intervals. 
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3 Z bosons 



These results allow for the establishment of the isolated photon sample as a baseline for further 
analyses, like for example, the production of photons accompanied by jets or, with larger data 
samples, constrain the nuclear PDF. 



3 Z bosons 



The Z bosons have been detected in CMS in their Z— > pi v\i and Z— >■ e + e decay channels. 
In the following, results in the muonic decay channel, based on an integrated luminosity £ = 
7.2 pb _1 of PbPb collisions data, are presented. The selection in this channel exploits their quite 
clear signature: two oppositely charged high transverse momenta muons (pj > 10 GeV/ c). 
In this way, 39 Z candidates are extracted, with an almost negligible background. Figure [3] 
presents the dimuon invariant mass. The Z mass peak in PbPb data has a resolution similar to 
that obtained in pp data. 
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Figure 3: Dimuon invariant mass spectra in PbPb collisions. Full squares are opposite-sign 
dimuons, while the empty circle shows a unique like-sign dimuon candidate. The histogram 
shows the corresponding distribution measured in pp collisions at 7 TeV within 60-120 GeV/ c 2 , 
scaled to the 39 PbPb candidates. 



Several kinematic distributions have been studied with this sample of Z candidates, like the Z 
rapidity and px, shown in Figure[4]left and right, respectively. They are compared to Powheg |T3l 
predictions (perturbative QCD at NLO) for pp — > Z — > scaled with the number of binary 

collisions, and also to other models simulating PbPb interactions which incorporate nuclear ef- 
fects such as isospin, energy loss and shadowing |T4l Il5l The experimental distributions are 
consistent with them all, not being still precise enough to quantitaively distinguish between 
them. 

When attending to the centrality of the PbPb collision, no dependence of the production of Z 
bosons on the number of participants is observed, as shown in Figure [5] Using the prediction 
for pp — > Z — > from Powheg, scaled with the number of binary collisions, and the 

minimum bias yield of Z — > in PbPb interactions, a relation Raa = 1-00 ± 0.16 ± 0.14 is 

established, validating the Z bosons as a reference particle unaffected by final-state interactions 
in the strongly dense system produced in Pb-Pb collisions at the LHC. . 
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Figure 4: Yields of Z — > fi]d per event in PbPb collisions: Left: dN/ dy versus Z boson rapidity, 
y. Right: d 2 N/ dydpx versus the Z boson pj. Vertical lines (bands) correspond to statistical 
(systematic) uncertainties. See References section for the different nuclear PDF models. 
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divided by the expected nuclear overlap function Taa- Vertical lines (bands) correspond to sta- 
tistical (systematic) uncertainties. See References section for the different nuclear PDF models. 
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4 W bosons 



4 W bosons 

W bosons have been detected and measured at the LHC for the first time in heavy ion collisions, 
in their W ± — > p^v decay channel. Their distinctive experimental feature is the presence of a 
single high px muon recoiling in the transverse plane to a significant imbalance in the total 
momentum of the event, which accounts for the undetected neutrino. 

Figure [6] (left) shows the muon transverse momentum distribution (solid dots) for events se- 
lected online with a low (2-3 GeV/ c) muon p T threshold and fulfilling some quality cuts, such 
as a minimum number of hits used in their offline reconstruction or requiring a distance in the 
transverse plane between the muon impact parameter and the primary vertex to be less than 
300 }im. The muon pseudorapidity region studied is |?/^| < 2.1. A Z-veto is applied based on 
the presence of a second high-pj muon and the resulting dimuon invariant mass. An enhance- 
ment of events is observed in the region p T > 25 GeV/ c, where the contribution from muons 
from W and Z decays is expected. In fact, a fit (black solid line) to the data is performed, 
taking into account two contributions, one for the W— > ]iv signal (modelled with a Pythia pp 
— > W — > }iv simulation, the normalization left as a free parameter in the fit) (green-hatched 
histogram) and another for the background, described by a 3 parameter function (blue-dashed 
line). The fit describes the data quite well. 

The imbalance in the transverse momentum in the event, {tfj), is computed as the opposite 
sign of the vectorial sum of all charged particle transverse momenta in the event, with pj > 3 
GeV/ c. The mean value of this quantity for all online selected muon events is shown in Figure|6] 
(right) (black squares) as a function of the event centrality. When requiring the presence of a 
high-px muon (p T > 25 GeV/ c) the associated mean value of j/j shifts to higher values of 
~ 40 GeV/ c, and is far less dependent on the centrality of the collision. This result agrees 
with that expected for the typical tfj produced by an undetected neutrino originating from W 
decay. To enhance the contribution from the W signal, events are therefore required to have 
p T > 25 GeV/c and ffj > 20 GeV/c. 

With this selection, the transverse mass distribution, given by mj = J 2p T j/x(l — cos <p), being 
(p the difference in azimuthal angle between the muon and tfj vectors, of the remaining events 
is presented in Figure [7] for PbPb (red dots). The data sample of pp collisions collected at 
the same centre-of-mass energy is analyzed in a similar way, and their mj distribution is also 
shown in the figure (blue-open squares). They both represent compatible W signals, almost 
background free, as they are in good agreement with a W Pythia signal simulation, embedded 
in PbPb events generated with Hydjet (PYTHIA + HYDJET [16]) (green-hatched histogram). 
Residual contamination from other electroweak processes (Z— > and W ± — > t ± v) (2%) 

are subtracted, and QCD background is estimated to be ~ 1% and included as a systematic 
uncertainty. 

From the sample of W candidates selected, 275(301) W+ and 264(165) W+ in PbPb(pp) data, 
the W production cross section in PbPb collisions, scaled by the nuclear factor Taa and in pp 
interactions, per unit of W rapidity, is computed and shown in Figure [8] (left) as a function of 
the number of participants in the collision, N par t- Two main effects are observed: one, the W 
production in PbPb does not depend on the centrality of the collision and the other, the com- 
parison of the individual W + and W~ production in PbPb and pp data, reflect the different u 
and d quark content in Pb nuclei and in the proton. The Raa factors obtained are: Raa(W + ) = 
0.82 ± 0.07 ± 0.09 and Raa(W~) = 1.46 ± 0.14 ± 0.16. Nevertheless, the total W (W+ plus W~) 
production is consistent, at LO, with N co u scaling, as Raa(W) = 1.04 ± 0.07 ± 0.12 

The difference in W + (ud — > W + ) and W~ (dti — > W~) production at LHC and the spin con- 
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Figure 6: Left: Single-muon transverse-momentum spectrum for \rj^\ < 2.1 in PbPb data (red 
points). Signal (green-hatched histogram) and background (blue-dashed histogram) contribu- 
tions are fitted (black-solid line) to the data. Right: Mean value of tfj for charged tracks as a 
function of the event centrality, before any selection on muon-triggered data (black squares) 
and after it (red circles), together with predictions from PYTHIA + HYDJET [16J (green trian- 
gles). 
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Figure 7: Transverse mass distribution for selected events in PbPb (red-filled circles) and pp 
(blue open squares) data, compared to PYTHIA + FFYDJET simulation (hatched histogram). 
The error bars represent statistical uncertainties. 
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References 



servation rule in their leptonic decay originate a different angular yield of W + — > }i + v and 
W~ — >■ }i~v, known as charge asymmetry and given by A = (%+ — Nw-)/ (Nw+ + Nw-), 
being Nw the efficiency corrected yield of W candidates. The distribution of the asymmetry as 
a function of the muon is presented in Figure [8] (right) for both PbPb and pp collisions at 
= 2.76 TeV, together with predictions from MCFM HZ| generator using MSTW08 HO PDF 
and in addition, for the PbPb data, the EPS09 HU nuclear PDF. 
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Figure 8: Left: Centrality dependence of normalised W — > yiv cross sections in PbPb collisions, 
for all W candidates (red filled points) and separated by charge, W + (violet-filled squares) and 
W~ (green-filled stars). The open symbols at N par t = 120 represent the minimum bias events. 
At Npart = 2, the corresponding cross sections are displayed for pp collisions for the same 
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Right: Charge asymmetry as a function of muon 1 1] \ for PbPb (red-filled circles) and pp (open 
blue squares) collisions. Overlaid are predictions from MCFM calculations with MSTW08 PDF 
(blue curve) and additionally with EPS09 PDF for PbPb collisions. 



As a conclusion a new measurement has been provided of the W production cross sections at 
a different centre-of-mass energy in the pp colliding system, validating once more the Stan- 
dard Model prediction, as observed in Figure [9] (left) where the W production cross section is 
depicted as a function of y/s. For the first time, the W yield has been measured in PbPb in- 
teractions, establishing the W as a reference particle unaffected by the hot and dense medium 
created, as seen in Figure [9] (right), where the Raa nuclear factor for the W as for the rest of 
electroweak bosons (Z and isolated photons) are shown to be consistent with unity over a wide 
range of energies, establishing these particles as "clean" probes of the initial state of the colli- 
sion. 

Whenever the uncertainties (both statistical and systematic) will allow it, the analysis of these 
bosons will be very useful to constrain predicted nuclear effects (shadowing, energy loss) and 
give access to quark/ antiquark/ gluon parton density functions (PDF) in the proton and neu- 
tron. 
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